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We have computed the contribution of zero modes to the value of the 
number of particles in the model of discrete (2+l)-dimensional nonlinear 
Schrodinger equation. It is shown for the first time that in the region of 
small values of the Chern-Simons coefficient k there exists a universal 
attraction between field configurations. For k = 2 this phenomenon 
may be a dynamic origin of the semion pairing in high temperature 
superconducting state of planar systems. 
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The cooperative behavior based on universal topological features of planar 
systems is the subject of a number of papers |], HI- The topological properties 
of (2 + l)-dimensional systems are described by a Chern-Simons term in 
the Lagrangian of the model. This term displays the phenomenon of the 
chiral invariance violation in such systems. The magnetic Chern-Simons 
field leads usually || f|, |], H 1] to the effective repulsion between field 
configurations. This property was described in detail in Refs. [[5|, |TU| where 



the structure of zero modes in continuous || and spacial discrete |T0[ models 
of (2 + 1)D nonlinear Schrodinger equation was studied. The competition of 
the basic nonlinearity corresponding to attraction, diffraction and additional 
nonlinearity, describing the repulsion caused by Chern-Simons interaction, 
resulted in the increase of the critical value of the number of particles N in 
the region of small values of coefficient k. The increase of N represents the 
existence of additional repulsion. 
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We argue in this paper that on the (2 + 1)D lattice there is universal at- 
traction due to Chern-Simons correlations between field configurations and 
show the conditions when this phenomenon takes place. We shall show that 
consideration of complete contribution of the statistical Chern-Simons fields 
in the form of the holonomies causes the additional attraction (to the bare 
one) at small numbers k of links. The necessary condition for that is nat- 
urally the existence of bare attraction between field configurations. This 
result is the first indication that there exists the universal attraction due to 
Chern-Simons correlation. The comparison with the results of our previous 



papers |9|, |10[ shows that the condition for obtaining the attraction is the 
inclusion in the consideration the Wilson and Polyakov exponents on the 
equal footing. In other words, we consider the model which takes into ac- 
count the conditions of the compactness for the temporary component of the 
gauge potential as well as for its spatial component. The gauge invariance 
requirements lead immediately that discrete evolution should be considered 
(for detail see jnj). We want to emphasize that the gauged discrete (2 + l)D 
nonlinear Schrodinger equation gives us a convenient tool to display this 
(2 + 1)D system phenomenon of additional attraction which has a general 
and universal character. 

The equation of motion in the model of the gauged discrete (2 + 1)D 
nonlinear Schrodinger equation has the form 

(4 + ty + h.C. - 4)pm,n = ~1C (? m n - p m) n $m{w m , n - 1) . (1) 

Here t x = e tA m, n +d/dm j g operator of so-called magnetic translations. 
The parameter C = g \k\ in Eq.(l) contains the coupling constant g of 
the classical nonlinear Schrodinger equation and Chern-Simons coefficient k. 
Besides the consideration of the model on 2D lattice [|l(J we included the 
discrete time t 6 2 in the description. This leads for the stationary states 
^(m,n,t) = p m . n e %t to the existence of nonlinearity presented by the sine 
function in Eq.(l). The multiconnection of the 2D manifold has been taken 
into account by the gauge field A^(m, n) = {w m>n , A^ n , A^^) where 

W m>n = [( A 2 G ( m -m',n- OX/4'X + Pm',n'+l) A m',h'- 

m' ,n' 

(AiG(m -m',n- n')){p 2 m , , + P^+inOA^n'] (2) 
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is the temporal component of Chern-Simons potential and 

Arh,n = A 2 G(m -m',n- n')p 2 m , n , (3) 

m',n' 

is the x-component of the vector potential. The notation Am,n denotes that 
the components A x (m,n) are determined on the links connecting the sites 
(m,n), (m+l,n). In Eqs.(2),(3) Ai >2 /(r) = /(r + ei j2 ) — /(r) is the gradient 
on the lattice with coordinates r = {m, n) G Z 2 ; is the unit vector. In 
accordance with the rules of the gauge field theory on the lattice we assumed 
that while the gauge field is defined on the lattice links, the curl of the field 
A M (r) and the density p 2 are defined on the sites of the dual lattice. Green 
function on the lattice in Eqs.(2),(3) has the form: 

/" d 2 k e *{M m - m ')+ fc i/( n -™')} _ i 

G(m — m'n — n')= / - — — — . (4) 

v ' 1 J (2vr) 2 4-2cosfc :c -2cosA; y v 1 

The main goal of this paper is to study the dependence of the critical 
number of particles N = J2 Pm'.n' on the parameter C = g\k\ considering the 

ml ,n' 

arbitrary large contribution of the temporal component of the gauge potential 
as well as of its spatial component. To solve the problem we compute this 
dependence using the zero modes found for various values of Chern-Simons 
coefficient k. 

We performed the simulation of the problem (l)-(4) on the lattice with 
linear size L < 20 using the method of the stabilizing Petviashvili multiplier 



12|| . The description of the block diagram of the method and the details 
of calculations may be found in Refs. |], |10j- The form of the functions 
p m ,n, Arh,n, u! m ,n found numericaly is displayed in Ref . [pToH . The form of these 
functions in the present paper is qualitatively the same. We used the function 
p mtn for the calculation of the dependence N(C) which is shown in Fig. 1. 
Note that the limit k — > oo is equivalent to the zero contribution of the gauge 
fields A m n , ?/; min , when in the continuous limit N cr = 11.703. 

This number of particles separates the 2D collapse regime at iV > N cr 
and its absence at iV < N cr . On the lattice the values Nftj? are always 



less than the critical number of particles in the continuous limit even if 
Chern-Simons fields are neglected. Therefore the main problem of interest is 
whether N^ t (A IJj ^ 0) is smaller or greater than N l c f(A^ = 0) if the Chern- 
Simons gauge fields are taken into account. 
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Figure 1: The dependence of the critical number of the particles on the 
parameter C for the three types of the nonlinearity in Eq.(l) (see for details 
the text). The dashed line shows the value N lat (A^ = 0) = 11.605. 
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Figure 2: The lines of equal value of the field p m n for C = 1 (a) and C = 5 
(b). 

The curve "a" in Fig.l shows the result of the paper [10] when we consid- 
ered only the part 4 — 2 cos A x — 2 cos A y of the contribution to Eq. (1) of the 
spatial gauge field components and did not consider the discrete time. The 
case "b" in Fig.l corresponds to complete consideration of spatial gauge field 
contribution to the l.h.s of Eq.(l) with the same properties of the time as it 
is above. The dependence N(C) in the case "c" in Fig.(l) presents the result 
of computations on the (2 + 1)D lattice with consideration of the complete 
contribution of all (arbitrary large) gauge field components taking into 
account the discrete time. 

From dependence N(C) in the case "c" in Fig.(l) one can infer that 
the decrease of the critical number of particles N with decreasing of the 
parameter C is equivalent to the increase of the attraction in comparison 
with the case when we do not take into account the contribution of the 
Chern-Simons fields. 

It is known |], Q that in Chern-Simons systems there arises the induced 
angular momentum proportional to 1/k. The calculated lines of the equal 
value of the field p m ^ n for C — 1 are shown in Fig. 2. Here we would like 
to pay attention to the extended s-symmetry of the field p mi n for a small 
value of the parameter C = 1 at great distance from the origin. We found 
weak display of this phenomenon (see Fig. 2). The accuracy of the performed 
calculations according to our estimates is several percents. The comparison of 
the results on Fig. 2 for the different value of the parameter C shows that the 
display of extended s-wave symmetry increases with decrease of parameter 
C. 
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The origin of the considered phenomena is as follows. We consider the 
nonlinear nonlocal dependence of the components A^m, n) of Chern-Simons 
gauge field via the field p m>n in its complete form presented in Eq.(l) for 
the Laplacian. In particular, if we extract (JTOj a part of this contribution 
to nonlinearity, specifically taking into account the compact version p(4 — 
2cosAj. — 2 cos A y ) of the nonlinearity p(A 2 x + A y ) 0, the rest part in the 

discrete Laplacian in the continuous limit has the form cos A x -^+cos A y -^ ■ 
The decrease of the coefficients in this expression in comparison with the 
unity for the value of A m<n in the region of small coefficients k leads to the 
decrease of diffraction. It is seen that the origin of additional attraction due 
to Chern-Simons fields is the decrease of the diffraction. Neglecting this effect 
we observe |10[ only the repulsion due to Chern-Simons fields. The anisotropy 
of this operator is the reason of extended s-symmetry of the ground state. 

As it was pointed out above the considered phenomenon exists on the 
lattice under the condition of the complete taking into account the gauge 
field by the holonomies e* Afl . The physical interpretation is clear: because 
the zero component of gauge potential plays the role of chemical potential 
the arbitrary large value of the one corresponds to arbitrary large value of 
the energy which we add to the system adding a particle. 

Note that discrete models are characterized by the features which are 
absent in the continuous limit. In particular, the localized states can exist 



in (1 + 1)D discrete nonlinear Schrodinger equation [[L3|l . In our case the 
dimensionality of the problem as well as the discrete character of the space 
and the time are important. Using the arguments in inverse order we have to 
include in the model discrete space and time on the equal footing in order to 
consider the large magnitude of the gauge field in the form of Wilson exponent 
as well as Polyakov exponent having in mind the gauge invariance. Our 
simulation shows for example that without the condition of the discretization 
of the time we can not obtain the attraction due to Chern-Simons fields. 

Finally, we should like to make a general remark. The attraction be- 
tween particles in the systems with Chern-Simons interaction was a subject 
of extensive studies during the last ten years. The attention was focused 
[fL4l , [15 , [16(1 on analyzing the symmetry of the state (s-, p- or d-wave states) 
with non-zero value of the superconducting gap in the framework of the per- 
turbation theory when the parameter a = 4n/\k\ was small, i.e. in the limit 
\k\ ^> 1. The discovery []I7| of the time-reversal symmetry-breaking p-wave 
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superconductivity in S^RuO^ (see also Refs. [|T8|, |XQ| , p0| ) stimulated the re- 
cent paper |H] where the search for the induced Chern-Simons term in P- 
and T-violating superconductors was performed. 

The picture of the Chern-Simons correlations, described in the present 
paper, is in some sense beyond the above-mentioned approaches. We found 
the attraction due to the Chern-Simons gauge field in the essentially non- 
perturbative region of the small value of the coefficient k irrespectively of the 
symmetry of the ground state. The choice of the coefficient k itself determines 
the symmetry of the considered state. In this respect, Fig. 2 presents only one 
of the possible symmetry of the ground state characterized by the specific 
value of the parameter C. 

Let us suppose in the Eq.(l) that C = 0. We obtain in this case the 
model describing the nonlocal interaction of Chern-Simons vortices. This 
case corresponds in continuous limit to universal nonlinearity of p 5 kind (from 
the point of view of scale transformations) and differs from Gross-Pitaevskii 
model f%| when the nonlinearity in the equation of motion is a local one 
and proportional to p 3 . We plan to study this interesting limit in a separate 
paper. 

In conclusion, we studied the dependence of the critical particle number 
on the link numbers of the field configurations. Using the model of the 
discrete (2 + 1)D nonlinear Schrodinger equation we found for the first time 
the existence of the attraction due to Chern-Simons fields. It was shown that 
the origin of this phenomenon is the suppression of the free propagation by 
Chern-Simons fields at small link numbers. Note that for g = 1 the semion 
value k = 2, which is of topical interest, lies inside this region. Therefore 
the found attraction may be the dynamical reason of the semion pairing and 
phase transition to a superconducting state. 

We are grateful to S.N. Vlasov, E.A. Kuznetsov, A.G. Litvak, A.M. Sa- 
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